Ballistic transport in ferromagnet— superconductor— ferromagnet trilayers with 

arbitrary orientation of magnetizations 
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Transport phenomena in clean ferromagnet-superconductor-ferromagnet (FSF) trilayers are stud- 
ied theoretically for a general case of arbitrary orientation of in-plane magnetizations and interface 
transparencies. Generalized expressions for scattering probabilities are derived and the differential 
conductance is computed using solutions of the Bogoliubov-de Gennes equation. We focus on size 
and coherence effects that characterize ballistic transport, in particular on the subgap transmission 
and geometrical oscillations of the conductance. We find a monotonic dependence of conductance 
spectra and magnetoresistance on the angle of misorientation of magnetizations as their alignment is 
changed from parallel to antiparallel. Spin-triplet pair correlations in FSF heterostructures induced 
by non-collinearity of magnetizations are investigated by solving the Gor'kov equations in the clean 
limit. Unlike diffusive FSF junctions, where the triplet correlations have a long-range monotonic 
decay, we show that in clean ferromagnet-superconductor hybrids both singlet and triplet pair cor- 
relations induced in the F layers are oscillating and power-law decaying with the distance from the 
S-F interfaces. 

PACS numbers: 74.45.+C, 74.78.Fk 



I. INTRODUCTION 

The interplay between ferromagnetism and supercon- 
ductivity in thin film systems is a phenomenon that 
attracts considerable interest of researchers for some 
time already. 1 Apart from potential device applica- 
tions, ferromagnet-superconductor (F-S) hybrid struc- 
tures can be major tools in the quest for novel su- 
perconducting phenomena^^ Variety of interesting the- 
oretical predictions, such as the existence of 7r-state 
superconductivity in F-S multilayer systems^*^ and 
characteristic oscillations in the superconducting tran- 
sition temperature as a function of thickness of ferro- 
magnetic layers (2&2ii2iiiiiS have been already confirmed 
experimentallyjiiiiiiSiiSiiLi^ while the existence of the 
spin valve effecliiSiSSi was more difficult for confirmation 
than expected. 21,22 Other peculiar phenomena that have 
been predicted to exist in F-S hybrids with inhomoge- 
neous magnetization - the occurrence of a long-range 
triplet pairing(2i2iSSi2Si2i or the so-called inverse prox- 
imity effect^ - still wait for experimental realization. 
Hence, a proper understanding of these effects is a prereq- 
uisite for setting up appropriate experimental conditions. 

Long-range triplet pairing, in particular, was first pre- 
dicted as a consequence of proximity of an inhomoge- 
neous ferromagnet to a superconductori2ii2!i2£i2& How- 
ever, it was shown that spin-triplet pair correlations may 
also arise in layered structures consisting of supercon- 
ductors and homogeneous ferromagnets with differently 
oriented magnetizationsiS2i2S The simplest structure of 
this kind is an FSF system with homogeneous but non- 
collinear magnetizations of the ferromagnets. In the par- 
allel (P) or the antiparallel (AP) alignment such corre- 
lations are absent. It was found that in diffusive junc- 
tions triplet correlations have a long-range monotonically 
decaying component in the ferromagnetic layers^ Very 
recently, effects of these correlations seem to have been 



observed in half-metallic ferromagnets ,2L2ii22i 

In this paper we study a variety of effects that occur 
due to the interplay of ferromagnetism and superconduc- 
tivity in clean FSF double junctions. It has been already 
established for the case of collinear magnetizations that 
the most significant consequences of quantum interfer- 
ence in such structures are the subgap transmission and 
the periodic vanishing of the Andreev reflection^ The 
annulment of the Andreev reflection occurs at the ener- 
gies of geometrical resonances, which correspond to the 
maxima in the density of states. Here, we investigate 
these phenomena in FSF hybrids with an arbitrary angle 
between magnetizations. In particular, we address the 
issue of triplet correlations, and describe its influence on 
transport properties. 

We use solutions of the Bogoliubov-de Gennes equa- 
tion in the scattering formulation (Section II) to obtain 
the probabilities of processes that charge carriers undergo 
and calculate the conductance spectra for arbitrary ori- 
entation of magnetizations. Section III illustrates some 
typical cases. In particular, the conductance spectra of a 
clean FSF hybrid with a thin superconducting layer (i.e., 
such that the layer thickness d s is less or of the order of 
the superconducting coherence length £o = frv-p/nAo) 
are dominated by subgap transport since majority of 
the charge carriers is transferred directly from one elec- 
trode to another without interaction with superconduct- 
ing condensatei2ii2i On the other hand, situation in the 
junctions with a thick superconducting layer {d s £o) is 
quite the opposite: majority of the carriers form Cooper 
pairs, thereby increasing the subgap conductance; also, 
the conductance spectra show pronounced oscillatory be- 
havior above the gap due to interference-generated geo- 
metric resonances inside the superconductor. However, 
we show that no extraordinary effects arise when the rela- 
tive orientation of magnetizations is between parallel and 
antiparallel - the spectra for intermediary angles simply 
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fall in between those for P and AP alignment. The ob- 
tained results allow us to infer some general conclusions 
about the nature of ballistic transport properties of FSF 
hybrids, which we illustrate on the example of depen- 
dence of magnetoresistance on voltage and S-layer thick- 
ness (Section IV). In Section V we address the problem 
of spin-triplet pair correlations when magnetizations are 
non-collinear. By solving the Gor'kov equations in the 
clean limit, we show that both triplet and singlet super- 
conducting correlations induced in ferromagnets are os- 
cillatory functions of distance from S-F interfaces, with 
a power law decay. This result is consistent with ob- 
tained monotonicity of transport properties as the angle 
of relative orientation of magnetizations is varied. 



II. THE MODEL 

We consider an FSF double junction consisting of a 
clean superconducting (S) layer of thickness d s , con- 
nected to clean ferromagnetic layers (F) of thickness df 
via thin insulating interfaces (I), Fig. 2] To describe the 
ferromagnets we use the Stoner model for inhomogeneous 
exchange field h(r) lying in y-z plane, parallel to the lay- 
ers. The exchange field has, in general, different orien- 
tations in the left and right ferromagnet, which (without 
loss of generality) we take to be symmetric with respect 
to the z-axis 



h(r) = ho (0, ± sin(a/2), cos(a/2)) , 



(2.1) 



for | a; | > d s /2. The case a — (a = it) then corresponds 
to the P (AP) alignment of magnetizations. 
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FIG. 1: A schematic of an FSF junction with magnetiza- 
tion vectors lying in y-z planes, at angles ±a/2 with the z 
axis. The a:-axis is perpendicular to the layer interfaces. The 
strength of insulating barriers (I) at the interfaces is measured 
by a dimensionless parameter Z — 2Wo/hvY- 



Quasiparticle propagation is described by the 
Bogoliubov-de Gennes equation 



for the four-vectors \I>(r) = (wf( r ) u l( r ) v l( r ) v \( r )) > 
where the Hamiltonian can be written in the compact 
form as 



H = ? 3 ® (So (r) - h(r) • ff) + A(r) ® 1. 



(2.3) 



Here, % and Oi are the Pauli matrices in orbital and spin 
space, respectively, 1 is 2 x 2 unity matrix, 



#o(r) = 



ffo(r) 






#o*(r) 



(2.4) 



H (r) = -h 2 V 2 /2m + W(r) + U(r) - fi being 
the one-particle Hamiltonian, U(r) and fi arc the 
Hartree and the chemical potential, respectively, while 
E is the energy with respect to /i. The inter- 
face potential is modeled symmetrically by W(r) = 
Wo [S(x + d s /2) + 8{x — d s /2)], where S is the Dirac 
delta-function. The matrix A is given by 



A(r) 



A(r) 
A*(r) 



(2.5) 



where A(r) is the pair potential. We will not consider 
the case of a ferromagnetic superconductor(25iii£i2L2£ and 
henceforth we take h(r) = outside the ferromagnets 
and A(r) = outside the superconductor. Since there 
is only one S layer in the system, without loss of gener- 
ality we can choose A to be real. The electron effective 
mass m is assumed to be the same throughout the junc- 
tion. For simplicity, the Fermi energy of the supercon- 
ductor and the mean Fermi energy of the ferromagnets 
are assumed to be the same, E-p. The model could be 
straightforwardly generalized to the case of a Fermi ve- 
locity mismatch between the layers (see, e.g. Ref . l33l) . 

Wave-vector components have to be normalized by us- 
ing the condition for canonical transformation, which af- 
ter integration over the volume of the S layer reads 

J d 3 r[K(r)| 2 + K(r)| 2 ] =1. (2.6) 

Here, a denotes the spin orientation, (a =1,1) and a 
is opposite to a. As the exchange field removes spin 
degeneracy, a spin-generalized gap equation has to be 
used. 39 Using the condition for the singlet pairing in the 
superconductor, 

K (r>* (r) - u x (v)v] (r)] = 0, (2.7) 



for \x\ < d s /2, the gap equation can be written in familiar 
form 



A(r) = V(r)J2 M r K( r ) + u l( r K( r )] tanh 



E 



U9(r) = EV(t) 



(2.2) 



(2.8) 

In Eqs. (|2.7(l and (|2.8J) the summation is performed over 
wave vectors in the superconductor q, taking into ac- 
count the dispersion relation between q and E. Also, we 
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assume for simplicity that the pairing interaction poten- 
tial V(r) = V = const, inside the superconductor. 

The parallel component of the wave vector, kn ff) is 
conserved due to translational invariance of the junction 
in directions perpendicular to the x-axis, but differs for 
the two spin orientations. Consequently, the wave func- 
tion can be written in the form 



/u T (x)e 8K i 



*(r) = 



M|(a;)e ik li.i- 
u;(a;)e ik il.T- 
\ «f (x)e* k n-i' 



\ 
) 



(2.9) 



In the following, we will use the stepwise approximation 
for the pair potential, A(r) = A&(d s /2 — x)Q(d s /2 + xV 
where denotes the Heaviside step function. In Ref. |33 
the averaged self-consistent pair potential is calculated 
and can be used instead of the bulk value. Fully self- 
consistent numerical calculations have been performed 
recently for FSF 12 - 40 and NS 4 i geometries (here N stands 
for a normal non-magnetic metal). 

The scattering problem for Eq. I|2.2[) is solved assum- 
ing bulk ferromagnetic layers, df — ► oo. Then, the four 
independent solutions of Eq. 11'. 21 correspond to the four 
types of injection involving an electron or a hole incidence 
from either the left or the right electrode^ 

Let us first consider the electrons (E > 0) with positive 
x-components of the velocity. Solutions of Eq. (|2.2|) are: 
for the left F layer (x < —d s /2) 



u x (x) 



cos(a/2)e lk t x + cos(a/2)6 T e- < ^ a: 



-ism(a/2)b ie -* k t x , 



-t sm(a/2)& T e~' ifc T - + cos(a/2)e ik t 3 



- cos(a/2)6je~ 



-ik^ " 
ikf 2 



= —i sin(a/2)a|e zfc T x + cos(a/2)a-f e lk i x , 



cos(a/2)ajy 



i sin(a/2)a-f e 



for the S layer (|x| < d s /2) 



— / \ - iqT x , - — iqTx . — iq^ x . - — iq~ x 

Uf(x) = uc\e + uc2e H ^ + vc^e ^ + vc^e H ^ , 

— / \ - iq, x , - — iqi~ x , — iqTx . - —iqTx 

u l\ x ) = ucse^i + uc^e y ± + vc^e H ± + vc%e y ± , 



vi(x) = vcie lq i x 



+ uc 3 e tq i x + wc 4 e iq f 



— / \ - iqTx i - —iqTx . — iq . x , - —iq, x 

w-)-(x) = vc^e^i + vc^e *i +ucre^i + uc%e , 
and for the right F layer (x > d s /2) 



u T (x) 



cos(a/2)c^e k T a 



■ i sin(a/2)c|e lfc J- 1 



ui(x) = isin(a/2)cie ik r x +cos(a/2)cie ik t x , 



v r [x) 



zsin(a/2)d i e" lfc T x + cos(a/2)<i T e lk i x , 



cos(a/2)d i e~ lfc T x + ism(a/2)d^e i 



Here, u = v / {l + fl/E)/2 and v = - fl/E)/2 are 
the usual BCS coherence factors, and fl = \/E 2 — A 2 is 



the quasiparticle kinetic energy with respect to the Fermi 
level. Perpendicular (x-) components of the wave vectors 
in the F layers are 



fc± = ^{2m/t?) (E F + p a h a ±E)- k 2 ^, 
while in the S layer they are given by 



<Z± = J(2m/V) (E F ±Q)-kl <j , 



(2.10) 



(2.11) 



where p a = 1 (—1) for a =1 (j). The superscript (+ or 
— ) corresponds to the sign of quasiparticle energy. 

At the layer interfaces the wavefunction ^(r) is con- 
tinuous, while its first derivative has a discontinuity pro- 
portional to a dimcnsionless parameter Z = 2Wq/Hvf 
measuring the height of potential barriers at the inter- 
faces, 



*( r )U=±(<i s /2 T o) 
<9*(r) 

dx a:=±(d s /2=F0) 



*(r)U=±(d s /2±o), (2-12) 
<9*(r) 

dx x=±(d B /2±0) 

-k F ZV{r)\ x=±ds/2 . (2.13) 



This yields a system of 16 linear equations in as many 
unknowns. By solving it, we find the amplitudes a a , b a , 

Ccn dai Ci, . . . , Cg. 

The scattering of an electron coming from the left leads 
to four possible processes: (1) local Andreev reflection 
(i.e., formation of a Cooper pair by two electrons from the 
left F layer); (2) normal reflection; (3) direct transmission 
to the right F layer; (4) crossed, or nonlocal, Andreev 
reflections^ (i.e., formation of a Cooper pair by two 
electrons from the opposite F layers). The respective 
probabilities are given by 



Re 



tan 



B a = \b a 

C a = |Cr 

D„ = Re 



tan 2 - Re 4 !fe 



tan 



- Re ^ 



tair 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



Conservation of probability current gives the usual nor- 
malization condition 



A a + B a + C a + D 



1. 



(2.18) 



Solutions for the other three types of injection can be ob- 
tained by symmetry arguments. In particular, the calcu- 
lated probabilities should be regarded as even functions 
of EM 
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It can be shown that A. 

d 



D a — whenever 
{it -Q7)= 2n7T 



(2.19) 



for n = 0,±1,±2, . . ., independently of /i an d Z. There- 
fore, both direct and crossed Andreev reflection vanish 
at the energies of geometrical resonances in quasiparticle 
spectrum. The absence of direct and crossed Andreev 
processes means that all quasiparticles with energies sat- 
isfying Eq. I|2.19|l will pass unaffected from one electrode 
to another, without creation or annihilation of Cooper 
pairs. The effect is similar to the over-the-barrier res- 
onances in a simple problem of one-particle scattering 
against a step-function potential^ the superconducting 
gap playing the role of a finite-width barrier wL 

Both the presence of insulating barriers and exchange 
interaction reduce A a and C a and enhance B a and D a . 
Approaching the tunnel limit (Z — > oo), the spikes in 
A a , C a , and D a , as well as the dips in B a , occur at the 
energies given by the quantization conditions 



d s q+ = 7Ji7r, d s q a = n 2 Tr, 



(2.20) 



which correspond to the bound-state energies of an iso- 
lated superconducting film. In this case, these bound 
states are the only conducting channels, both for super- 
current and quasiparticle current 



III. DIFFERENTIAL CONDUCTANCES 



1.0 
0.8 



55 



0.2 



^0.6 
§0.4 
0.2 
0.0 



' 1 I 


l < l 

Z=0 














(a) 






Z=l 








(b) 



FIG. 2: Differential conductance spectra of an FSF junction 
with (a) Z = and (b) Z = 1, for a thin S layer (d s k F = 10 3 ), 
X — 0.5, and three relative orientations of magnetizations: 
a = 0, a = it/2, and a — n (top to bottom). 



When voltage V is applied to the junction symmetri- 
cally, the charge current density can be written in the 

fo rr ~ 34.48.49 



J(V) 



(2-Kfh 



x J A\ la (1 + A a -B a + C a -D a )8f(k,V) 

OO 

= _4_ f dE Y P CT x 

J z-f 

x J d 2 k|| iff (A CT + C7 ff )5/(k,y), (3.1) 

where P a = (1 + Pa X)/2 with X = h /E F , and Sf(k, V) 
is the asymmetric part of the nonequilibrium distribution 
function of current carriers. In the last equality the nor- 
malization condition, Eq. I|2.18|l . was taken into account. 
Without solving the suitable transport equation we take 
Sf(k, V) = f (E-eV/2)-f (E + eV/2), where /„ is the 
Fermi-Dirac distribution function^L^ In this approach, 
the charge current per orbital transverse channel is given 
by 



I(V) 



dE [f (E - eV/2) -f (E + eV/2)} G(E), 

(3.2) 



where 



G(E) = G ]T P a J (A a + C a ) 



(3.3) 



is differential charge conductance at zero temperature 
and Go = e 2 /irh is the conductance quantum. 

One can see that Eq. I|3.3|) is a simple generalization 
of the Landauer formula^ with terms that take into ac- 
count transmission of the current through Cooper pairs 
and quasiparticles. For E < A, the subgap transmission 
of quasiparticles (without conversion into Cooper pairs) 
suppresses the Andreev reflection, while for E > A all 
the probabilities oscillate with E and d s due to the in- 
terference of incoming and outgoing particles. 

The influence of the exchange interaction and rela- 
tive orientation of magnetizations on the conductance 
spectra is illustrated for X = 0.5 and Z = 0, for thin 
(d s k F = 10 3 , Fig. 0) and thick (d s k F = 10 4 , Fig. EJ} S 
films. In all the illustrations the bulk value of supercon- 
ducting pair potential is characterized by Aq/E f = 10~ 3 , 
which corresponds to ^o^f = (2/-7r)-10 3 . For thin S layers, 
d s k F < 10 3 , the bulk value Ao is replaced by the aver- 
aged self-consistent pair potential, calculated in Ref. |3^. 
The spin-polarized subgap transmission of quasiparticles, 
and consequently a strong suppression of the Andreev 
reflection, are significant in junctions with thin S films, 
whereas the conductance oscillations above the gap be- 
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FIG. 3: Differential conductance spectra of an FSF junction 
with (a) Z = and (b) Z = 1, for a thick S layer (d s fc F = 10 4 ), 
X = 0.5, and three relative orientations of magnetizations: 
a — 0, a — 7r/2, and a — tt (top to bottom). The spectra of 
the corresponding FIS junctions (generalized BTK results^) 
are shown for comparison (dashed lines). 



FIG. 4: Zero-bias conductance of an FSF junction with Z — 
(top panel) and Z — 1 (bottom panel) at zero temperature, as 
a function of the normalized S layer thickness d s /£o, for X = 
0.5, and three relative orientations of magnetizations: a = 0, 
a — 7r/2, and a = tt (top to bottom). The corresponding NSN 
values (X = 0) are shown for comparison (dashed lines). 



come pronounced as d s is increased. The magnetoresis- 
tance is apparent, as conductance is greater for the P 
(a = 0) than for the AP (a = tt) alignment. The pres- 
ence of non-collinear magnetizations (a = tt/2) does not 
lead to any nonmonotonicity in the conductance spec- 
tra with respect to collinear cases (a = and a = tt) - 
the values for a — tt/2 simply fall in between the P and 
AP curves. Moreover, it can be seen that the conduc- 
tances for E > A oscillate in phase ioi a — 0, a = tt/2, 
and a = tt, Fig. [Ja). For thin S films (d s k F < 10 3 ) 
transmission of the spin polarized current and suppres- 
sion of Andreev reflection are still dominant at energies 
below A, Fig. Ufa). The conductance spectra exhibit 
similar behavior for finite transparency of the interfaces. 
This is illustrated for weak non-transparency [Z = 1, 
Figs.|3b) andC3»]. From both Figs. |1 and it can be 
seen that conductances attain their high-energy values, 
corresponding to conductances of an FNF double junc- 
tion, when E is of the order of several A. 



IV. MAGNETORESISTANCE 

In clean FSF structures electrons propagating from F 
to S layer are not necessarily converted into Cooper pairs. 



Moreover, part of them has always a non-zero probabil- 
ity of being directly transmitted from one F electrode 
to the other, even at the voltage below 2A/ei^24 This 
probability increases with decreasing d s . As the applied 
voltage is increased, the conductance spectrum starts to 
resemble the one of an FNF junction. Fig. 0] shows the 
zero-bias voltage G(0) as a function of normalized thick- 
ness of the S layer, d s /£o, for three relative orientations 
of magnetic moments: a — 0, a — tt/2, and a = tt. 
It can be seen that mutual differences between the cor- 
responding values of G(0) are larger the thinner the S 
layers are. At thicknesses of the order of the supercon- 
ducting coherence length (<i s /£o ~ 1) these differences 
slowly disappear, and vanish completely at d s /£o ~ 5. 
The increase of the normal reflection probability with d s 
is more rapid than the increase of the Andreev reflection 
probability if the insulating barriers were present in the 
junction. Consequently, zero-bias conductances are more 
sensitive to d s in the junctions without (Z = 1) than in 
those with the full transparency (Z = 0), as can be seen 
in Fig. H 

For hybrids with thin S layers (d s k-p ~ 10 3 or less) 
zero-bias conductance decreases monotonously when a 
changes from to tt (Fig. [5}, both for Z = and Z = 1. 
If the S layers are thick (d s kp ~ 10 4 or greater), G(0) no 
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longer depends on a: the Andreev reflection dominates 
at energies close to the Fermi level, while the probability 
of direct transmission from any of the F layers through 
the superconductor is negligible, independently of mutual 
orientation of magnetizations. At energies above the gap, 
dependence of G(E) on a is also monotonous, with values 
close to those for the corresponding FNF junction. 




FIG. 5: Zero-bias conductance of an FSF junction with Z = 
and Z = 1 at zero temperature, as a function of the relative 
orientation of magnetizations a, for X = 0.5 and thin S layer, 
d s k F = 10 3 . 
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FIG. 6: Magnetoresistance of an FSF junction with (a) Z — 
and (b) Z = 1 at zero temperature, as a function of the 
applied voltage, for X = 0.5, and ultrathin, thin, and thick S 
film, d s k F = 10 2 , 10 3 , 10 4 , respectively. 



Magnetoresistance is defined as 

-Rap — Rp 



MR 



(4.1) 



where -Rp(AP) = V/-fp(AP) is the junction resistance in 
P (AP) alignment of magnetizations. In Fig. El MR is 
shown as a function of bias applied symmetrically at the 
ends of an FSF junction with Z — and Z = 1, for 
X = 0.5 and three values of d s k F : 10 2 , 10 3 , and 10 4 . 
The main characteristic of junctions with thin S layers 
(d s /cF ~ 10 2 — 10 3 ) is the dominance of direct transmis- 
sion over Andreev reflection. Also, the average number 
of quasiparticles converted to Cooper pairs is small, inde- 
pendently of voltage, making the number of up spins that 
cross from one F layer to the other much greater than the 
corresponding down spins. This difference is greater in 
the P than in the AP alignment of magnetizations, lead- 
ing to a magnetoresistance that gets more pronounced 
as the S layer becomes thinner. Conversely, in junctions 
with a thick S layer (d s kp ~ 10 4 ), direct transmission 
probability is very low and hence most of the electrons 
that enter the S layer at energies below A are converted 
into Cooper pairs of net spin zero. Thus, for the elec- 
trons in one F layer the influence of the opposite F layer 
becomes practically negligible, and therefore MR ~ for 
eV/2A < 1. At higher voltages, eF/2A > 1, MR rises 
practically monotonously due to a gradual increase of di- 
rect transmission. 

Therefore, magnetoresistance of clean FSF junctions 
exhibits a strong dependence on the S layer thickness, 
even at a low bias. The MR is always more pronounced 
if the exchange field is stronger. However, it decreases 
towards zero with increasing d s due to the dominance of 
supercurrent over the normal current. 



V. SPIN-TRIPLET CORRELATIONS 

A peculiar property of proximity of singlet-pairing su- 
perconductors and inhomogeneous ferromagnetic metals 
is inducement of triplet correlations between electron- 
hole pairs of equal spin orientations^ These correlations 
also exist in systems with locally homogeneous exchange 
fields of non-collinear configuration. 27 Triplet correla- 
tions in F-S hybrids only resemble those of magnetic su- 
perconductors: while in the latter the Cooper pairs could 
be truly spin-triplet, in the former they remain singlet de- 
spite the presence of both singlet and triplet components 
of anomalous Green's function. 

To study the triplet correlations in clean FSF junctions 
we solve the Gor'kov equations generalized to take into 
account non-collinearity of magnetizations in the F lay- 
ers. In this case, the net Green's function G(r 1 ,r 2 ) is a 
(2 X 2) ® (2 X 2) matrix in the Nambu space. The Gor'kov 
equations can be compactly written as 



(El - H) <7(ri,r 2 ) = 5{r x - r 2 )l, 



(5.1) 
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where Ti. is given by Eq. (|2.3(l . while i is the unity ma- 
trix in the Nambu space. As the system of our interest 
is part-by-part homogeneous, in each layer of the junc- 
tion the components of the matrix Green's function Q 
are functions of the relative coordinate r = ri r2- It 
consists of four matrix blocks, 



Q = 



G+ F 
F* G~ 



(5.2) 



Here, we are interested only in the anomalous block 



F = 



F ll F lt 



(5.3) 



The antisymmetric combination (F\± — F^)/2 describes 
singlet correlations, while the symmetric functions (-Ffl + 
F^)/2, Fff, and F±± are related to triplet correlations. 
If the exchange field h(r) has the same direction in the 
left and the right F layer (a = 0), then without loss 
of generality we can choose this to be direction of the 
z axis, see Fig. ^ The blocks Hq(t) — h(r) ■ a in H 
are then diagonal, and thus only the singlet correlations 
are present. The function F = F^i = F^ is then the 
usual pair amplitude, while F^ = F^ = 0. In general, 
however, it can be the case that < a < ir. Then, the 
functions Ffi and F^ are different, while and 
are non-zero, leading to existence of triplet correlations. 
For ferromagnets, where A = 0, Eq. 1)5. lfl reduces to 



lEl-H + (r))F(r)=0. 



(5.4) 



General solution of this equation for the right F layer are 



Fiji = cos(a/2)C 1 e lk t x + cos(a/2)C 2 e 



F u - 



+i sin{a/2)C 3 e lk t x + i sin(a/2)C 4 e- ife ^(,5.5) 
i sm(a/2)de^ x + i sm(a/2)C 2 e- lk t x 
+ cos(a/2)C 3 e lk t x + cos(a/2)C 4 e~^ x , (5.6) 
F iT = ism(a/2)C 5 e lk r x +isin{a/2)C 6 e- lk i x 

+ cos(a/2)C 7 e lk i x + cos(a/2)C & e-' lk i x , (5.7) 
F TT = cos(a/2)C 5 e lk i x + cos{a/2)C 6 e~ lk T x 

+tsm(a/2)C 7 e lk ^ x + % sm(a/2)C&e~ lk i x (£.8) 

modulo phase factor exp(ik|| CT • r). Solutions in the left 
F layer are obtained if we substitute a by —a, and use a 
different set of constants, Cg, . . . , C\§. To find the com- 
plete set of unknown constants C± , . . . , C\e we have to 
use the continuity of F aa i and its derivative at the S-F 
interfaces. However, to apply the boundary conditions 
at x = ±d s /2 we need to find solutions for F aa i inside 
the superconductor, where G acr ' and F aa t are coupled 
through the pair potential A(r). To remain within a 
tractable analytical procedure, we use a simpler method 
which is sufficient for qualitative discussion. First, we 



eliminate eight out of sixteen constants using condition 
that holds at the outer boundaries, 



F n 



x=±(d s /2+d f ) 



0. 



(5.9) 



Then, we connect these functions with corresponding so- 
lutions in the S layer, which for this purpose we treat as 
unknown parameters, defined as 



F n 


x=±d B /2 


= A U 


(5.10) 


F n 


x=±d B /2 


= A 2l 


(5.11) 


F^ 


x=±d B /2 


= 5i, 


(5.12) 


F ll 


x = ±d B /2 


= B 2 . 


(5.13) 



Note that the functions F aa > depend on E, kn, and x. In 
order to describe the spatial variation of triplet compo- 
nents we introduce the pair amplitudes 



f a<r ,{x)=N{Q) 



dE 



d 2 k 



||,<7 



27rfc| 



F aa -'(E,ku „,x) 



(5.14) 

where N(0) is density of states at the Fermi level. Using 
these functions, we can further construct 3 



fo(x) 
fs(x) 



[/U(3 



1 

2 



fu(x), 



(5.15) 



The following symmetries have to hold: fo(—x) = fo(x), 
fi(-x) = -fi(x), and / 3 (-a;) = f 3 (x). The function f 3 
describes singlet, while fo and /i describe triplet corre- 
lations in the system. The function fo corresponds to 
electron-hole pairs of zero net spin orientation, while f± 
corresponds to pairs of net spin orientation equal to ±1. 
The f 3 component exists even when the exchange field 
is absent. Both f 3 and fo fall off rapidly in a ferro- 
magnet: characteristic decay lengths are of the order of 
£f = ^JhDf I ho in diffusive, and £f = hv-p/ho in bal- 
listic heterostructures, where Df is a diffusion constant 
of a dirty ferromagnet. In diffusive FSF junctions these 
components are of a short range, since the exchange field 
tends to align spins, while f\ is of a long range and mono- 
tonically decaying. 52 Here, we will argue that in clean 
FSF hybrids f\ can be of the same range as f 3 , and 
no long-range monotonically decaying triplet correlations 
are generated. 

If we assume that the F layers are bulk, df ^S> d s , and 
focus on correlations in the vicinity of the S-F interfaces, 
\x\ > d s /2, then any plain wave of the form e ±lk " df will 
be rapidly oscillating with respect to e ±lk ° x . Therefore, 
the most important contribution to the integral of F al7 i 
over E in Eq. (|5.14|) is from the energies close to E — 
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0. The energy dependence in the wavevectors is thus 
negligible, and we can write k^ ~ kf(E = 0) = k a . 
Within this approximation the number of undetermined 
coefficients, defined by Eqs. <|5 . 1U|> — (|5 . 13|> . can be further 
reduced by applying the symmetries /j-f (jc) = fn(x) and 
fli(x) = —fii{x) that now hold. Hence, A\ = — A2 = A 
and B\ = B2 = B. Performing the integration over 
kn a in Eq. (|5.14|) . we finally obtain the following set of 
expressions for the pair amplitudes of the right F layer 

fo(x) = 0, (5.16) 
fx{x) = Bl(x) + isinaAJ(x), (5.17) 
f 3 (x) = Al(x)-ismaBJ(x). (5.18) 

where x = x — d s /2 is the relative distance from the right 
S-F interface. In Eqs. (|5.16l) - (|5.18[l we have introduced 
the auxiliary functions 



l(x) 
J(x) 



Qt(20 



Si (5 



Qi(x)-Qi(x), 



where 



Qa{x) 



d 2 k||^ sin [k<j (df — x)] 



2-Kkl 



sin (k a df) 



(5.19) 
(5.20) 

(5.21) 



The function l(x) is even and has a maximum at the in- 
terfaces, 1(0) = 1, while J(x) is odd, with J(0) = 0. 
The expressions for the pair amplitudes for the left 
F layer are straightforwardly obtained by substitutions 
a — ► — a and x — » — x. For a = tt/2 the auxiliary func- 
tions, Eqs. (|5.19(l and (|5.20J) . then simplify to 



l(x) 
J(x) 



sin(Ka:) 



sin(«;a;/2) 



(kx/2) 



sin(ftx) cos(kx) 
(kx) 2 kx 



(5.22) 
(5.23) 



where n = fcpf — k-p [ is the difference between the Fermi 
wavenumbers for the two spin subbands, 



(1 + X) 1 ' 2 -(l-X) 1 ' 2 } , 



(5.24) 



and, as before, X = ho/E-p. The approximated functions 
l(x) and J(x), given by Eqs. (|5.22|) and l|5.23|l . are shown 
in Fig. □ for X = 0.5. 

Under the assumptions made, the unknown parame- 
ters A and B in Eqs. (|5.17|l and (|5.18|l will depend on 
X and a. If X = and/or a = 0, ir then B = 0, since 
the triplet correlations are absent in these cases. From 
Eqs. (|5.1()|) - (|5.18|l we can draw the following conclusions 
concerning triplet correlations in clean FSF heterostruc- 
tures. The function /o is approximately zero inside the 
F layers. Since this function can only be generated by 
the exchange field it also cannot exist in the S layer. 
The singlet correlations, captured in f^, consist of two 
terms. If a = or 7r, fa reduces to Al(x). This term 




FIG. 7: The functions l(x) (solid curve) and J(x) (dashed 
curve), under approximation df 2> d 3 , shown for a — n/2 
and X = 0.5. 



has a decay length 27t/k ~ Thus, the range of singlet 
correlations increases with decreasing strength of the ex- 
change field. In special case of an NSN junction (X = 0) 
we obtain a well-known result that singlet pair ampli- 
tude monotonically decays into the N layerS For 1^0 
and non-collinear magnetizations two components of the 
triplet pair amplitude fi oscillate on the same scale and 
with the same decay length. 

We conclude that in clean FSF junctions at zero tem- 
perature both singlet and triplet pair correlations in fer- 
romagnets are oscillating and power-law decaying with 
the distance from S-F interfaces. Physical intuition be- 
hind this result is that in such systems there exists only 
one characteristic length determining decay of correla- 
tions in the ferromagnets, which weakly depends on ex- 
citation energy E. This stands in a contrast with find- 
ings for diffusive FSF hybrids where two characteristic 
lengths are present: and £ w = y / Df/2h\uj\, where 
uj = (2n + l)7rfcBT. At temperatures just below the 
critical one for the superconducting phase transition, 
T < T c , the two lengths set up two different scales, 
£uj ~ \/fiDf /2irkBT c ^> £/. Note that £ w is also the 
length that determines the decay of singlet component in 
nonmagnetic normal metal. 



VI. CONCLUSION 

We have studied the influence of misorientation of 
magnetizations on the properties of ballistic transport 
in clean FSF trilayers. By solving the Bogoliubov-de 
Gennes equation we have derived generalized expressions 
for probabilities of processes that charge carriers un- 
dergo. We use these probabilities to compute differential 
conductances for arbitrary orientation of magnetic mo- 
ments and interface transparencies. Preferability of di- 
rect quasiparticle transmission to the Andreev reflection 
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in thin S layers and more prominent resonant oscillations 
in thick S films are the main consequences of quantum 
interference in clean heterostructures. The subgap con- 
ductance is larger for P than for AP alignment as a result 
of strong magnetoresistive effect in thin S layers: when 
the S layer thickness is less or comparable to the super- 
conducting coherence length the direct transmission of 
spin polarized quasiparticles across the superconductor 
becomes a dominant transport mechanism. However, we 
show that no extraordinary effects arise when the rela- 
tive orientation of magnetizations is between parallel and 
antiparallcl - the spectra for intermediary angles simply 
fall in between those for P and AP alignment. 

From the obtained results we can further draw some 
general conclusions about the nature of ballistic trans- 
port in clean FSF trilayers with inhomogeneous magne- 
tizations. Firstly, the zero-bias conductance depends on 
the relative orientation of magnetizations only in het- 
erostructures with thin superconducting layers, decreas- 
ing monotonously when tuning from P to AP alignment. 
Secondly, magnetoresistance displays qualitatively differ- 
ent voltage dependance for thin and for thick supercon- 
ducting layers due to different quasiparticle spectra which 
tends to gapless or bulk, respectively. Thirdly, magne- 
toresistance increases with the thickness of the S layer 



and vanishes eventually. The effect is more pronounced 
for strong exchange fields and transparent interfaces. 

Non-collinearity of magnetizations leads to formation 
of spin-triplet pair correlations in FSF structures. Un- 
like the diffusive case, where the triplet correlations have 
a long-range and monotonically decaying component, we 
have shown that in clean ferromagnct-superconductor hy- 
brids these correlations are oscillating and power-law de- 
caying with the distance from S-F interfaces, similarly to 
the usual singlet correlations. This similarity in behavior 
of singlet and triplet pair correlations induced in ferro- 
magnets is the main reason why the transport properties 
of clean FSF junctions have monotonic dependence on 
the angle between magnetizations. These findings sug- 
gest that any spectacular features of triplet correlations 
are highly unlikely to occur in clean FSF structures. 
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